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We provide a microphysical model that connects neutrino oscillations to dark energy, which has 
predictions of Lorentz and CPT violating neutrino oscillations. We argue that the DE is a BCS 
condensate of flavored neutrinos. As neutrinos propigate in their own condensate they naturally 
have oscillations proportional to the DE energy density. All that is assumed in this model is a 
covariant coupling of neutrinos to gravity and a finite number density of neutrinos in the cosmic 
rest frame; this situation yields an attractive channel for the formation of a spin zero neutrino 
condensate leading to late time acceleration self consistently. Moreover, the vacuum oscillation has 
two predictions: (1) CPT violating oscillations proportional to the DE density (2) A connection 
between the evolution of the equation of state of the DE condensate and the neutrino oscillations. 
These predictions can be probed independently with future Supernovae and Neutrino Telescopes. 

PACS numbers: 



I. INTRODUCTION 

While we know observationally that dark energy (DE) 
makes the universe accelerate on large scales and might 
require a modification of General Relativity in the 
IR, nothing is known about how it may interact with 
matter. Rectently Ando, Kamionkowski and Mocioiu 
(AKM) proposed a new way to account for DE-induced 
Lorentz/ CPT violation — introducing neutrino oscilla- 
tions that are energy independent [5] • Moreover, ever 
since the discovery of DE a number of coincidences have 
been identified 1 . A peculiar coincidence exists between 
the observed neutrino mass differences between propi- 
gation states Am^ tm ~ 10 _3 eV 2 and the DE density 
Pde ~ (lCP 3 eV^) 4 . Some authors have exploited this fact 
by seeking to make connections between the evolution of 
the neutrino density and the quintessence scalar, which 
leads to promising phenomenology [H H]. However, in all 
models that connects neutrino physics to dark energy, a 
quesitons still remains: 

• Does neutrino oscillations and the emergence of 
DE come from the same underlying physical phe- 
nomenon? 

In this paper we attempt to answer the above question by 
providing a mechanism that non-trivialy ties in neutrino 
oscillations with the emergence of dark energy. The idea 
is straightforward: Due to a generic four fermion coupling 
from general relativity a flavored condensate of neutrinos 
form a scalar bound state (Cooper-pair) in the early uni- 
verse. At late cosmological times, the coupled equations 
of the neutrinos and the scale factor generically drives 
late time acceleration. Likewise, the effective action of 
the condensate interacting with Dirac neutrinos predict 
two types of neutrino oscillations: 

• Mass oscillations without the need of a see saw 
mechanism 1 , that arises from interaction between 



neutrinos and the dark energy condensate. 
• Lorentz and CPT violating oscillations that are in- 
dependent of energy in contrast to a E~ 1 decay for 
mass oscillations. 
The paper is organized as follows: In section II we derive 
the four-fermi interaction from general relativity coupled 
to free fermions followed by the evaluation of the effective 
potential of the neutrino condensate. Section III gives a 
description of the mechanism of flavor oscillations. We 
conclude with a discussion of precision electroweak tests 
of the model and future research. 



II. FIELD EQUATIONS 

Our starting point is to discuss how a gravitationally in- 
duced neutrino BCS condensate forms in the early uni- 
verse; due a weak attractive four-fermion interaction and 
a Fermi surface. We assume that the early universe is 
filled with a finite number density of neutrinos and follow 
the dynamical consequences. We will find that conditions 
for condensation arise naturally in general relativity. In 
fact general covariance plays an important role in the 
emergence of the four fermion interaction. A free fermion 
field subject to a covariant coupling with the gravita- 
tional field induces non-propagating torsion. When we 
solve for torsion field we naturally get a four-fermion 
contact interaction induced in the Einstein-Hilbert ac- 
tion. Due to the symmetries of FRW space-time a self- 
consistent Cooper-pairing between fermions of opposite 
spin is naturally induced. 

To address exactly how the four-fermion interaction 
arises we must introduce two independent fields to grav- 
ity: the tetrad, e„, an orthonormal coordinate basis for 
each point on the manifold, and a spin connection ui^ 



1 This mechanism assumes no Majorana mass terms and a van- 



ishing coupling constants of the Higgs to Dirac neutrino mass 
terms. 
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which connects (parallel transports) the tangent spaces 
at different points of the manifold. Note that lower-case 
Latin letters starting from fi,v,... denote spacetime in- 
dices, while capital Latin letters starting from I, J, . . . 
denote internal Lorentz indices. One can then associate 
a 4-dimcnsional spacetime metric g^ v via 
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(11.1) 



where the Minkowski metric is to be viewed as the met- 
ric of the internal space. Internal indices are raised and 
lowered with the Minkowski metric, while spacetime in- 
dices are raised and lowered with the spacetime metric. 
The requirement that the spin connection be torsion free 
is simply wj^j = 0. Let us now rewrite the Einstein- 
Hilbcrt action in terms of the spin connection and the 
vierbcin. 

The Einstein-Hilbert action is given by the well-known 
expression: 



Sfi — 
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d i xdet(e) 
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R(iv,e) - - (i^'e^V^ + cc.) 

(11.2) 

Note that a tetrad based formalism is essential for the 
inclusion of fermions in the theory, since Dirac spinors 
live naturally in SU(2). Therefore, covariant deriva- 
tives associated with the Dirac action are not the usual 
50(3, 1) covariant derivatives, but instead are given by 
V^ip :— d^j — jioj/'fijj ip, where -0 is a Dirac spinor, 
and j 1 are 4x4 gamma matrices. 

In the absence of fermions, the field equations for the 
connection gives rise to the metric compatibility (or zero 
torsion) condition: 
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which gives rise to the Christoffel symbols and the 
Einstein Hilbcrt action. However, things dramatically 
change when one includes the covariant coupling of free 
fermions to GR in the presence of first order GR. The 
interaction between the fermions and the connection 
changes the metric compatibility condition and the equa- 
tion for the connection. However, we can solve for the 
connection and metric compatibility by decomposing the 
connection into symmetric and antisymmetric pieces: 
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(II.4) 



,sIJ 



is the torsion-free spin connection satisfying 



when 

the compatibility condition (II. 3 1 and C 7J is the so called 



"contorsion" tensor. The idea is to integrate out the the 
contorsion tensor which then will lead us to the more fa- 
miliar second order formulation of gravity where connec- 
tions are just the metric dependent Christoffel symbols. 
This can be achieved simply by imposing the structure 
equations obtained by varying the action with respect to 
the connection. Using the identity 
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we can express the contorsion tensor in terms of the axial 

07 5 7 7 ?/), 



fermion current J5 



e^ JK = 4ttG 
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(II.6) 

From the above expression for the contorsion tensor it is 
clear that C m jK is a non-propagating field, its field equa- 
tions do not have any derivatives on it. Thus "integrating 
it out" is not only equivalent to reinserting its expression 



(III. 22 1 in the full action classically, but also quantum 
mechanically. Thus the four-fermion contact interaction 
term that we are going to generate in going from the 
first to second order formalism is quantum mechanically 
an exact result. This is a key difference from the "effec- 
tive" contact interaction that one obtains in non-abelian 
gauge theory where the mediating gauge fields do indeed 
propagate, and therefore the contact term is only a low 
energy approximation. 

Substituting (III. 22 1 into Hoist's action we find that the 
action can be written as [? 1 



S = S e [lj] + S d [uj] + Si, 



(11.7) 



The first and the second terms are the standard Einstein- 
Hilbert and Dirac actions involving Christoffel connec- 
tions. Crucially however, one obtains a third interaction 
term given by: 



Sint — 



d A x e J 5 j Jg 
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Such four-fermion interactions were already observed 
in Einstein-Cartan theory (■j 2 — > 00 limit in ( II.8[ >) , al- 
though they are suppressed by a power of Newton's con- 
stant (a factor of 1/k here). To form a condensate of 
fermions the Planck suppression must be transcended, 
which happens when the universe contracts to Planck 
densities, as we shall discover in the next section. 



A. Neutrino Superfluiditiy 

The key point in the last section was to understand how 
covariance leads to the four-fermion interaction when 
fermions interact with gravity. Various authors have real- 
ized that a weak and attractive interaction leads to a cos- 
mological BCS-like condensation of the fermions [T^] ■ 
The two known regimes where a four-fermion channel ex- 
ists are in the electroweak sector and in gravity. However, 
the electroweak case is repulsive for neutrinos 2 . On the 
other hand the gravitational case does gcnerically allow 



2 For a counterargument using modern arguments in fermi-liquid 
theory see 1 141 



3 



an attractive interaction and condensation is possible. In 
what follows we sketch how one obtains the "gap equa- 
tion." The four fermion interaction can be replaced with 
an auxiliary field, rendering the path integration as a 
gaussian over the fermionic field. This allows one to ex- 
plicitly evaluate the path integral involving the fermion 
self-interactions, leading to an effective potential for the 
condensate. The first step is to self-consistently find the 
generating functional for the condition that the auxiliary 
field A = (ip^ip). After expressing the spinors in the Weyl 
basis and introducing the auxiliary field, the four-fermion 
generating functional becomes: 



Z= [VA][D$[V{]e 



i(Sf or + Str. 
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(II.9) 

where we have integrated the Grassmann fields, defined 
the effective action S c ff (often referred to as T in quantum 
field theory literature) and approximated the functional 
integral by the saddle point (mean-field approximation, 
[TO]). The effective action S c g can be evaluated by per- 
forming the Gaussian integrals in terms of fermionic co- 
ordinates. As usual, one ends up with a fermionic de- 
terminant which allows us to express the effective action 
(see e.g. QUEL]): 



S e ff — Strcc i 



d p 



ln(det Ap) , 



(11.10) 



where the tree action is obtained from the interaction 
action 

"W) 2 ~ 
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The determinant of A p is straightforward to compute: 



det A p = [a 



|p|+ M ) 2 -A 2 ][u 



(|p|- M ) 2 -A 2 ], 
(11.12) 

where A is the auxiliary field at the saddle point. 
The chemical potential is responsible for aligning the neu- 
trinos in the cosmic rest frame and is the underlying rea- 
son for the Lorentz violation. Accordingly, we are left to 
compute the effective potential 



Kff = -C 



M 2 



efif 



(11.13) 



where 



d 3 p 

(2tt) 3 



V(|p|+ M )2 + A2 + v /(|p|-M) 2 + A2 . 

(11.14) 



Plugging I into V c s above leads to 
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The above potential has a minimum given by the gap 
equation 



dA 
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The potential energy also includes a contribution from 
the chemical potential as well. We can obtain the chem- 
ical potential from noting that the total number n of 
fermions is given by |10j 



no = I d x etp"/ ip = 

■2 



Sn 
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(N+ 1 +lnA z 



a 3 n. 



(11.17) 



When the system relaxes at the minimum of its potential, 
the total gap energy density of the fluid is given by 
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The gap field A is a time dependent function and red- 
shifts as ~ &~ 3 , however at late times it tends to a con- 
stant, which is why it exhibits a late time de Sitter state. 



Solving the EH action ( II. 2 ) in an FRW background leads 



to the modified Friedman equations 
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where p m is an extra matter term and 
A 4 



Pgap 
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-(2N + 3 + 2hiA 2 ) . 



(11.19) 
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The cosmological evolution is represented as a transcen- 
dental equation that requires numerical treatment. The 
authors [H] were able to find that the modified Fried- 
mann equation leads to a late time attractor exhibiting 

acceleration (i.e. e = —jp < 1) that depends on the 
magnitude of the gap. In particular, when the gap has a 
mass M ~ 10~ 3 the gap is also mill-electron-volt scale, 
Aq ~ 10 -3 , which corresponds to an energy density with 
negative equation of state w ~ — 1. This condition leads 
to late time acceleration. In particular we can solve for 
the equation of state from the Raychaudhuri equation: 



_ lnp(A(t)) gap 

w(t)cff - hTo(^ ' 



(11.21) 



which is plotted in fig [T]. Hence, we have established 
that the neutrino condensate exhibits late time accelera- 



tion. 
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FIG. 1: This is a plot of the evolution of the dark energy 
equation of state parameter, u)(t) of the neutrino condensate. 
We see an early phase of Phantom evolution followed by a 
late time asymptote to a de Sitter phase, ie. u = — 1. 



III. THE MECHANISM 

We would like to explore how the condensate that results 
from the four fcrmion self-interaction induces oscillations. 
The key point is that the specific form of the four fermion 
operator in the presence of flavors naturally mixes the 
flavor indices due to the Fierz identity. When the neutri- 
nos of different flavor condensed we are lead naturally to 
flavor oscillations. First we need to extend the fermion 
interactions to (at least 2) flavors. It is possible to solve 
(II. 3 1 for N-flavors of fermions: 
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where B is a flavor index. This leads to a modified inter- 
action term, 



Ant = J d4x e ■0a757 I '0a'0b757i' ( /'b (III. 23) 
Flavors will get reshuffled due to a Fierz transformation: 



To preserve isotropy and homogeneity, the gap must 
evolve in a frame aligned with the the cosmic ex- 
pansion, hence the four- vector A M = (A, 0,0,0). dc 
Gouvea encountered a similar Lagrangian and inter- 
preted as a "matter potential", because it resembles 
electron neutrinos propagating in the presence of elec- 
trons where A is analogous to the electron num- 
ber density as neutrinos pass through the sun, i.e. 
(e^a^e) = {\f2GpN e , 0, 0, 0), where iV e is the electron 
number density [8 j . However our case represents different 
physics since the neutrinos are moving in the background 
of their own BCS-like condensate. From ?? we get the 
Dirac equations for a neutrino created from the conden- 
sate: 

~ \AB&»Wf a VaB ~ m\v% = , 



m A VaB = 0, (111.26) 



where v is a left handed neutrino field, a, a — 1,2 are 
spinor indices and A, B are the flavor indices. The above 
Dirac equation has, up to an extra 75, the same form 
as the CPT violating case considered by [5]. We now 
proceed to study the flavor oscillations in the presence of 
the condensate. 



IV. NEUTRINO OSCILLATIONS 

We would like to understand how the presence of the neu- 
trino condensate affects Majorana neutrino oscillations. 
For simplicity we will consider a two-flavor oscillation. 
The key idea is that the neutrinos of differing flavors can 
condense. A neutrino moving in this flavored conden- 
sate will experience an oscillation since the condensate 
has non vanishing off-diagonal components in the flavor 
eigenbasis. This analysis dovetails the MSW effect where 
the equation of motion of neutrinos in a background elec- 
tron density, N e , are taken into account. Recall that 
a neutrino in a flavor eigenstate can be decomposed in 
terms of its mass eigenstates \v-\) and \v 2 ). For concrete- 
ness consider an electron neutrino \v e ) 



\v e ) = cos9\v\) + sin6\v2) , 



(IV. 27) 



= V»A757 V'aV'b757i' ( /'b = ^aV'b^aV'b 

+i>Alh1pBi>Alhll>B + i'Al^B^Alli'B (IH-24) 



combining (II. 8 ) with (III. 24 1 and identifying the Dirac 
spinor with a neutrino field ip = v 

C = ivAd^^A + Pa7 m ^bA aiA b + (1 + 7 5 )A A b^a^b- 

(111.25) 

The second term is a Lorentz/CPT violating interaction 
between the DE condensate and the neutrino. The third 
term represents the DE-condensate acting like a Parity 
and Lorentz invariant masses, (1 + 7s)A = iha respec- 
tively. 



where 9 is the mixing angle that parameterizes the mixing 
matrix U. In general a flavor eigenstate is related to the 
mass eigenstate by 



where 



\v)a = U ai \v)i , 



TT ( cosO sinO\ 
M ~ \-sin9 cosO ) 



(IV. 28) 



(IV. 29) 



If a neutrino propagates as a plane wave, then its time 
development is 

v(t, x) e = cose e - ipi » x " 'i>i + sin9e- W2 ^% 2 ■ (IV.30) 
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An ultrarelativistic neutrino of mass m, traveling a dis- 
tance L along the z-direction will have a phase factor 
p^x^ = Et — px ~ (E — p z )L. We can simplify matters 



by realizing that, 
leading to: 

v(t, x) e = cosOe 1 



E 



P : 



E+p 



2E 



and E 



\P\; 



Now we are ready to describe the time development for 
neutrinos propagating in the condensate. 
We can transform from the mass into the flavor basis by 
making use of the that fact that Wlf = 1 and multi- 
plying both sides of ?? by Usi- In the ultrarelativistic 
regime, where the neutrino is traveling a distance L and 
the condensate is purely time-like, the equation of motion 
(III. 26 1 becomes: 



dt \v : 



Am 2 A ( 



sin 2 9 cos9sin0\ 
2E U \cos6sin0 cos 2 8 J 

(A + V2G F N e ) A ex /2 
AL/2 



(IV.32) 



Note the the electron number density N e is included to 
take into account for the MSW effect. v x is a linear com- 

\ ee — A xx . The equa- 



bination of & v T , and A = A ee — L 
tion for antineutrinos is identical to Eq. (IV.32) with 
Aab — > —Aab and N e — > — N e . One important differ- 
ence between the Lorentz preserving mass term and the 
Lorentz violating coupling is a factor E^ 1 in the mass 
coupling. 

This was pointed out in |2j and can have observational 
consequences at ultra high energies. In the absence of 
matter effects the oscillation probability is 



P ex = Pxe sin 2 26> eff sin 2 



A cff 



L 



(IV. 33) 



where 



A cff = y (Asin26 + A ex ) 2 + (A cos 20 - A) 2 , (IV.34) 
A cff cos 28 cS = A cos 26 - A, (IV.35) 
A off sin 26> cff = A sin 29 + A ex , (IV.36) 

and A = Am\/2E U . Henceforth A ex is assumed to 
be real. We immediately see that even when we set all 
masses m to zero, there will still be oscillations in terms 
of the dark energy condensate 



*cff 



(A ee - A„„) 2 . 



(IV. 37) 



In particular, the equation of state u> depends on the 
logarithm of the gap's density 



1 + u(t) e s 



]np(A(t)) 



gap 



In a(t) 



(IV.38) 



so we can obtain a correlation between the equation of 
state evolution uj c ff and the oscillation angle. 



Recently AKM performed an analysis of CPT violating 
oscillations [3 . These oscillations were parametrized by a 
four vector a p which has the same algebraic form as our 
condensate A M . Therefore we can map their parameters 
for Lorentz violation in to ours and use it to constrain 
the mass term obtained from the DE interacion rriA = A. 
We find that: 



m e f fsin2Q,ie lri = \ab-A 
m e f fsin2&de^ — Aab A* 



(IV. 39) 



where Xab is the coupling constant matrix, defined by 

x 



Xab A = A A b- AKM found that m eff ~ ^a°. Where 



SX is difference between the eigenvalues of the coupling 
matrix. In terms of our neutrino condensate the effec- 
tive mass is, m e // ~ 5AA. AKM found that future ex- 
periments that target cosmogenic ultra-high-energy (~ 
1CP 17 — !0~ 19 GeV) neutrinos produced by the interac- 
tion of ultra- high-energy cosmic rays with CMB photons, 
can improve the sensitivity to down to 6XA ~ 10~ 30 GeV. 
If we assume that ~ O(l) then the dark energy neu- 
trino condensate amplitude can be constrained from fu- 
ture high-energy neutrino observations. 



V. DISCUSSION 

In this note we have considered the effects of general rel- 
ativity covariantly coupled to neutrino fields. In an FRW 
background a gravitational contact attraction drives free 
neutrinos at finite number density into a BCS-like con- 
densate. Such a condensate has been shown to exhibit 
the observed late time acceleration provided that the con- 
densate has a milli-electron-volt VEV. In this work we 
showed that at late times if the condensate has off di- 
agonal flavor components neutrinos that propagate on 
cosmological scales will exhibit flavor oscillations, poten- 
tially connecting the scale of neutrino oscillations with 
the scale of dark energy today. 

One immediate concern is whether or not the repulsive 
interaction that exists between neutrinos and Z-bosons 
can prevent the gravitationally induced Cooper-pairing 
of neutrinos from forming. Likewise, if the condensate 
forms in the early universe, competing repulsive inter- 
actions and scattering can break apart the pairs. The 
condensate has correlations on the order of the Hubble 
radius at the time that it was formed. During the elec- 
troweak epoch tEw ~ 10~ 12 s, the repulsive neutrino self- 
interaction is mediated by the Z° boson. The effective 
4-point operator is: 



2M| 



V A a^VBV A a P v B 



(V.40) 



The weak bare coupling is clearly larger than the gravita- 
tional coupling. However, the condensation mechanism 
requires a coupling between neutrinos of different flavors. 
The Z-boson only couples neutrinos of the same flavor so 
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it will not compete with the gravitational attraction of 
neutrinos. 

Furthermore the the temperature of the early universe 
during the electroweak epoch is Tew ~ l0 3 GeV. There- 
fore, the SU (2) vector bosons are in equilibrium with the 
universe during this time. As a result the interactions due 
to W and Z bosons will be screened by the mechanism of 
Debye screening. The interaction potential for a SU(N) 
gauge theory with Nf fermions was found to be|15j: 



the condensate is j 



U(r)su(N) 



ran 



47rr 

r 



t)V*G f T, 



(V.41) 



Clearly during the electroweak epoch U(r, Tew)su(N) ~ * 
so the four-fermion electroweak exchange is suppressed 
in the early universe. The Debye screening will not af- 
fect the gravitational interaction since gravity is out of 
equilibrium at temperatures below the Planck scale. 
Another issue to investigate in the future is the stability 
of the condensate 3 . This issue is important since simi- 
lar MaVaN models that couple quintessence to neutrino 
gcncrically exhibit a linear instability with an imaginary 
speed of sound. The onset of the instability is close to 
the time that neutrinos become non-relativistic. One ex- 
pects that since our DE condensate couples to neutrinos 
in a similar manner to MaVaNs that an instability is in- 
evitable. 

The speed of sound c 2 s for an adiabatic perturbation of 

3 The issue of instability of the condensate is currently being in- 
vestigated by the author and Niayesh Afshordi. 
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3H(l + u) 



(V.42) 



In a future work we will analytically prove that pertur- 
bations of the fluid is stable (ie. the speed of sound is 
positive). The speed of sound is always positive (V.42 1 



U) is 



as we can see from fig [T] that the time derivative of 
positive and the denominator in the second term is neg- 
ative (since cu < —1). Therefore, hurestically we see that 
the speed of sound for perturbations of the condensate is 
positive and stable. 

There are a few things left for future investigation. First, 
we would like to understand how this type of oscillation 
might affect other sectors of the standard model such as 
flavor changing neutral currents. In particular, neutrino 
scattering with the condensate can enhance flavor oscil- 
lations. It will be interesting to evaluate these effects 
explicitly for future neutrino beam experiments. 
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